Nodal dynamics, not degree distributions, determine the 
structural controllability of complex networks 
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Structural controllability has been proposed as an analytical framework for making predictions 
regarding the control of complex networks across myriad disciplines in the physical and life sciences 
(Liu et al., Nature:473(7346):167-173, 2011). Although the integration of control theory and network 
analysis is important, we argue that the application of the structural controllability framework to 
most if not all real-world networks leads to the conclusion that a single control input, applied to the 
power dominating set, is all that is needed for structural controllability. This result is consistent 
with the well-known fact that controllability and its dual observability are generic properties of 
systems. We argue that more important than issues of structural controllability are the questions 
of whether a system is almost uncontrollable, whether it is almost unobservable, and whether it 
possesses almost pole-zero cancellations. 



INTRODUCTION 

How can we control complex networks of dynamical 
systems [THH]? Is it sufficient to control a few nodes, or 
are inputs needed at a large fraction of the nodes in the 
network? Which nodes need to be controlled? A recent 
paper [10] suggests that we can address these problems 
using the concept of structural controllability and 
in doing so we may be able to forge new connections 
between control theory and complex networks. The two 
main results from this analysis are (1) that the number 
of driver nodes, Nrj, necessary to control a network is 
determined by the network's degree distribution and (2) 
that Np tends to comprise a substantial fraction of the 
nodes in inhomogeneous networks such as the real- world 
examples considered therein. 

However, both conclusions hinge on a critical assump- 
tion of the model in [10]: the results (implicitly) re- 
quire that the "default" structures of the dynamical sys- 
tems at the nodes of the network have infinite time 
constants. This modeling assumption implies that, un- 
less otherwise specified by a self-link in the network, 
a node's state never changes absent influence from in- 
bound connections. However, the real networks consid- 
ered in [TO] — including food webs, power grids, electronic 
circuits, regulatory networks, and neuronal networks — 
typically manifest more general dynamics at each node, 
i.e. they typically have finite time constants [T2Tll4j . 

With this assumption, the minimum number of in- 
dependent control inputs required to ensure a technical 
property known as structural controllability |11) can be 
calculated for the network, as described in [JJj] . The main 
problem with the argument set forth in |10j is not a tech- 
nical one: indeed the assumptions therein are clear and 
the mathematical results are correct. Then, why are the 



results tenuous? Critically, structural controllability [TT] 
is premised on the idea that if the nonzero parameters 
in the mathematical model can be selected so that the 
system is controllable (an elementary concept in control 
theory; see for example |15j). then the system will be 
controllable for all parameters except a set of zero mea- 
sure. That is, if the system is controllable for one set 
of (initially nonzero) parameters, then controllability is 
guaranteed generically for that system. The results pre- 
sented in jlOj require that a critical assumption be made 
before applying the structural controllability approach. 
Specifically, it is assumed that each node has an infinite 
time constant. As we shall see in the next section, the 
assumption of an infinite time constant implies that a cer- 
tain parameter in the mathematical model of the system 
is equal to zero, and therefore that term is off-limits as 
far as structural controllability is concerned. As one can 
imagine, any approach to system analysis that only al- 
lows the modification of nonzero terms, makes the results 
potentially quite sensitive to which terms are set to zero 
in the first place. Indeed, if the infinitc-time-constant as- 
sumption is relaxed, and generic linear dynamics are as- 
cribed to each node, one obtains a categorically different 
result. Indeed, we show in this paper that all networks 
with finite-dimensional linear dynamics (save a special 
set of parameters of zero measure) are controllable with 
a single input. While mathematically true, such a con- 
clusion is neither reasonable nor practical for real-world 
networks, and thus calls into question the general ap- 
proach of applying structural controllability in this way. 

Assuming arbitrary (up to a set of measure zero) linear 
dynamics, we show here that (1) a single time-dependent 
input is all that is needed for structural controllability, 
and (2) that this input should be applied to the power 
dominating set (PDS) [TB] of the network. Thus for many 
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if not all naturally occurring network systems, structural 
controllability does not depend on degree distribution 
and can always be conferred with a single control input. 



MODELING NETWORKS FOR CONTROL 

Large interconnected systems are commonly repre- 
sented as complex networks [T7J [TB] . For many biological 
and physical networks, each node in the network corre- 
sponds to a dynamical system. Often, the dynamics of 
these nodes can be modeled by a system of ordinary dif- 
ferential equations [El [20] : 



N 

fc=i 



X>i«i(*), (i) 



where Xi is a state at node i, N is the number of nodes, 
P is the number of inputs, and the n 2 elements pop- 
ulate the adjacency matrix. Here, the term —piXi repre- 
sents the intrinsic dynamics at the node, absent external 
influences. The external inputs, Uj(t), enter the system 
through the coupling matrix {bij}. For analyzing control- 
lability, it is reasonable as a first step to consider purely 
linear dynamics as shown in Eq. ([I]) — an approach clearly 
articulated and well motivated by [TO] . 

Note that Eq. includes two terms in dynamics for 
Xi, one related to the linearization of the intrinsic nodal 
dynamics, namely —p^Xi, and one related to a potential 
self link in the model, namely anXi, related to the net- 
work topology. Although both terms are identical math- 
ematically, they arise from categorically different sources, 
and thus are not interchangeable. 

The term — pi is the pole of the linear dynamical sys- 
tem at each node, and Ti — 1/pi is the associated time 
constant. Rewriting in terms of transfer functions, we 
have 



Xi(s) = Gi(s) 



N P 
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where Xi(s) and Uj (s) are the Laplace transforms of state 
Xi(t) and input Uj(t) respectively, and 

Gi(a) = 



s + Pi 



is the transfer function of node i. This formulation is 
useful because it suggests inclusion of more general linear 
dynamics: the transfer function, Gi(s), can be replaced 
by any transfer function, of arbitrary order. 

The dynamics proposed in [TO] (see the supplemental 
material therein) are identical to ([2]), except that pi = 
for all i, namely Gi(s) — l/(s + 0) — a pure integrator. 
Written this way the simplifying assumption of the model 
in [TU] becomes clear: all subsystems by default have an 



infinite time constants (that is, the term pt — 0) un- 
less such dynamics are explicitly included in the network 
data set through nonzero diagonal terms, an ^ 0, in the 
adjacency matrix. 

However, infinite time constants at each node do not 
generally reflect the dynamics of the physical and biolog- 
ical systems in Table 1 of [10] . Reproduction and mor- 
tality schedules imply species-specific time constants in 
trophic networks. Molecular products spontaneously de- 
grade at different rates in protein interaction networks 
and gene regulatory networks. Absent synaptic input, 
neuronal activity returns to baseline at cell-specific rates. 
Indeed, most if not all systems in physics, biology, chem- 
istry, ecology, and engineering will have a linearization 
with a finite time constant. Thus while the model in [TU] 
does not proscribe self-links, this approach does place 
the onus on the modeler to ensure that any network rep- 
resentation includes such self-links where appropriate to 
compensate for the omission of the intrinsic nodal dy- 
namics that arise due to physical, biological, or other 
processes that, generally speaking, have nothing to do 
with network topology. 

To see the consequences of including generic nodal dy- 
namics on a network's structural controllability, we first 
rewrite the network dynamics in |2]) in state space form: 



■kit) = Ax(t) + Bu(t), 

A= [A-di&g{p 1 ,p 2 ,P3,-- 



(3) 
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where A € M ArxA ' is the adjacency matrix, and B E 
R NxP is the input matrix. The vector x(t) € R N is the 
vector of node states, and u(t) £ K p is the input vector. 

The system in Eq. |3|) is controllable if and only if the 
matrix 



[B, AB, ••• A N - 1 B] 



(4) 



is full rank, a standard result in control theory |15) . 
The system is said to be structurally controllable if the 
nonzero weights in A and B can be adjusted such that 
the matrix in Eq. Q is full rank [TT] . 

In [10] . the minimum number of driver nodes, Nn, 
is defined as the minimum number of inputs — i.e., inde- 
pendent, user defined, time-varying functions — such that 
when injected into the network guarantee structural con- 
trollability. This formulation explicitly allows each inde- 
pendent input to be connected to multiple (and possibly 
all) nodes in the network [101 ETJ . 

The paper [10] solves this minimum input problem us- 
ing an application of graph-theoretic concepts; their basic 
approach is to identify the number of "unmatched nodes" 
after finding a so-called maximum matching of the graph. 
Details are provided in the supplemental material of [TO] ; 
note also the prior analysis wherein the maximum match- 
ing theorem seems first to have been proved [35|. We 
observe that one can recast the poles at — pi as (nonzero) 
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self-links. But the set of all self-links [i — > i) is itself a 
maximum matching; all nodes in the network are then 
matched nodes. This implies that the network can be 
controlled with a single input, i.e. Nr> = 1, which follows 
directly from the maximum matching proof in [101 I23| . 

STRUCTURAL CONTROLLABILITY OF 
NETWORKS WITH GENERAL LINEAR 
DYNAMICS 

The following proposition provides a simple non-graph- 
theoretic proof that a "control hub" — a single driver 
node attached to all nodes — guarantees structural con- 
trollability with a single input. 

Proposition 1. For any directed network with nodal dy- 
namics in Eq. ^ (or equivalently Eq. [3]], with pi =/= 
and/or an ^ 0, i = 1, . . . , N, then Njj = 1. 

Proof. Select B = [1,1,..., 1] T (that is, connect a single 
input to all nodes). Lin's structural controllability theo- 
rem [TJJ states that if the system is controllable for one 
choice of the nonzero system parameters, then it will be 
controllable for all parameters except a set of measure 
zero. So, we explicitly construct a parameter set that 
makes the system controllable. Keep B as all ones, and 
choose Pi,P2i •• • iPn to be nonzero and distinct. Zero out 
all the network edges (i.e. nullify the adjacency matrix, 
A = 0). The system matrix A is now a diagonal ma- 
trix with distinct eigenvalues. Controllability of (A, B) 
follows by inspection. Thus, the system is structurally 
controllable and Njj = 1. □ 

By contrast the paper [TU] reported that for real- world 
networks, the minimum number of driver nodes Nd is 
strongly influenced by the sparseness and homogeneity 
of the network, as measured by the degree distribution, 
P(fcj n ,/c out ) (see [TU] for more details). Why did [TU] ar- 
rive at such different conclusions? Critically, the applica- 
tion of structural controllability does not consider vari- 
ations in system parameters that are a priori zero |llj . 
So, for example, if a link i j is absent, then ay = 0. 
The original paper [TU] allows for self-links but by de- 
fault does not include them. Further, the framework set 
forth in [10] assumes pi = (infinite time constant), and 
the network datasets in Table 1 of [10] do not include 
self-links to correct for this. Therefore, upon inclusion 
of first-order self dynamics, essentially all real networks 
are structurally controllable with Njj — 1, irrespective of 
network topology. 

In the case that the network topology does not explic- 
itly contain self links, the consequence of ascribing pure 
integrator dynamics (jpi = 0) to each node is categorical: 
the system is necessarily unstable. This is because the 
sum of the eigenvalues is given by the trace of the sys- 
tem matrix, which, in this case, would be trace(A) = 0, 



since there are zeros on the entire diagonal. This would 
imply that it is impossible to have any stable eigenvalues 
(negative real parts) without also having unstable ones 
(positive real parts), so that their sum is zero. There- 
fore, such a network of integrators must be purely oscil- 
latory or unstable, and cannot be asymptotically stable. 
Therefore, assuming pure integrators at each node, and 
no explicit self-links in the adjacency matrix, precludes 
passive stability which many natural systems enjoy. 

Have we taken the point about generic nodal dynam- 
ics too far? It may be desirable to model and control a 
network on a timescale that is faster than the dynamics 
of the intrinsic nodal dynamics. We concede that in such 
cases, it may be reasonable to treat the nodal dynamics 
as pure integrators (systems with infinite time constants) . 
However, we argue that structural controllability may not 
be appropriate for addressing these nuanced modeling is- 
sues. An essential feature of structural controllability is 
that no importance is assigned to specific values for the 
non-zero terms in the dynamics. Values are treated as 
either zero or not zero; there is no in-between. Thus, the 
choice of whether to zero out the self-loop terms a priori 
is a subtle modeling issue that should take into account 
the emergent timescales of the entire network. Therefore, 
we contend that model reduction [24] — which is essential 
for controller design — should be treated at the level of the 
entire network dynamics rather that at the level of indi- 
vidual nodes: indeed the timescales relevant for control 
are an emergent property of the system dynamics, and 
not strictly a feature of one node or another. With this in 
mind, we find that the tool of structural controllability — 
which is premised on a notion of generic parameters — is 
best suited to generic modeling assumptions. In this case 
this means assuming p t ^ 0, i — 1, . . . N . 

Above, we argue that structural controllability of com- 
plex networks depends on the dynamics at each node, and 
that only a single time varying input is required. Two 
questions remain: (1) How sensitive is structural control- 
lability to the dimension of the state space for each node? 
(2) Where should we inject the No independent time in- 
puts into the network, i.e. what is the minimum num- 
ber of nodes of the network to which the input must be 
connected? Proposition [I] explicitly depends on treating 
first order nodal dynamics as "self loops" in the network. 
Below we offer a more general treatment for arbitrary 
(linear) nodal dynamics that addresses both questions 
above. See Figure [IJ 

Given a directed graph, a PDS is, by definition, the 
smallest set of nodes such that all other nodes are down- 
stream of at least one node in the PDS. Obviously, con- 
trollability requires connecting the input (s) at least to 
this set; below we show that structural controllability is 
generically achieved by connecting a single input to the 
PDS. Before doing this, we need one definition: 

Definition 1. Suppose that there are K nodes in the 
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FIG. 1. Given a network, the PDS (large white circles) is the 
smallest set of nodes such that all other nodes (smaller grey 
circles) are downstream of them. Any network, with arbi- 
trary (and possibly different) order finite-dimensional linear 
dynamics at each node is structurally controllable from a sin- 
gle driver node (black square) tied to the PDS as shown. See 
Proposition [2] The edges in the structural control network 
are part of a minimum spanning tree (black edges, although 
this choice of edges, and indeed the PDS, is not necessarily 
unique). 



PDS. Attach a single control input, u, to this set via a 
control node. Augment the graph with this control node 
and add the K edges that connect it to the PDS. Then, 
all nodes are downstream of the input u (i.e. the control 
node is now the PDS of the augmented graph). Define 
the structural control network as an acyclic directed 
graph given by a directed spanning tree that starts at u 
and visits all nodes. 

We now state the main result. 

Proposition 2. Consider the nodal dynamics in 
with Gi(s) an arbitrary, proper, rational transfer func- 
tion 11 51/ of the form 



Gi(s) 



ni(s) 
di(s) ' 



where, 7i,(s) and di(s) are assumed to be generic polyno- 
mials ( all coefficients up to the order of the polynomial 
are assumed to be nonzero ) of finite but arbitrary order 
in s. Then, the network is structurally controllable with 
one (Nd = 1 ) independent input, connected to the PDS. 

Proof. Using the structural controllability argument, we 
are free to modify any nonzero parameters; if the sys- 
tem is controllable for one set of parameters, it will be 
gcncrically controllable. 

So, zero out all edges that are not in the structural 
control network and set all those in the structural con- 
trol network to 1; if this process results in a controllable 
system, as we now show it does, then the system will be 
controllable generically. 

All nodes in this structural control network are still 
downstream of u, but now there are no cycles. Since the 
structural control network is a minimum spanning tree, 
there is exactly one path between u and any specific node, 
i. Let J~i denote the set of nodes along the path from u to 



node i in the structural control network. Then transfer 
function from u to any given node is simply the product 
of the transfer functions along the path from u to the 
node: 



Hj(s) 



Xj(s) 
U(s) 



= J] G k (s), j = 1,2,. 



.N. 



(5) 



Since we may freely adjust the polynomial coefficients 
in the denominator terms, we do so to ensure there are no 
repeated poles in the entire network and similarly adjust 
the numerator coefficients to ensure no pole-zero cancel- 
lations along any path in the structural control network. 
Since there are no pole-zero cancellations, and all poles in 
the network are unique, a minimal realization of the N x I 
transfer function [Hi(s), i?2(s), . . . Hn(s)] t must contain 
exactly one eigenvalue for each pole of the network. It 
is obvious that the minimal realization requires no more 
eigenvalues than that. The number of eigenvalues in the 
minimal realization is equivalent to the number of eigen- 
values that are both controllable and observable. Thus 
all states are controllable for this parameter set and, by 
the' structural controllability theorem [IT], the network 
is structurally controllable. □ 

For first-order nodal dynamics, our main result is not 
substantively different from those presented for discrete 
time and finite state systems in [35] . They show that 
networks with nontrivial nodal dynamics are structurally 
observable with a single output node and structurally 
controllable with a single input node. Our modest gen- 
eralization to arbitrary-order nodal dynamics is at best 
incremental over their work. Indeed, the main contri- 
bution of our paper lies not so much in any technical 
advance as it does in providing a timely clarification of 

ma- 



SIMPLE EXAMPLE: A FOOD WEB 

To illustrate the ideas of this paper, consider a simple 
food web comprising one predator and one prey species. 
Let H denote the number of herbivores (prey) and C de- 
note the carnivores (predators). We begin by noting that 
historically, the classic models of predator-prey dynam- 
ics K7I take the form 



C = -jC + eCH 
H = aH — f3HC 



C = C{- 1 + eH) 
H = H(a-/3C), 



(6) 



where a, f3, 7, e > 0. Linearizing these dynamics about 
the nontrivial equilibrium (C* ,H*) = (^,7) this model 
has the following local dynamics: 



x = 







< V (. 

tL 
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where x = (SC, SH) is the vector of small displacements 
relative to the equilibrium (C* , H* ) . Note that the lin- 
earized food web is fully connected (whereas [10] include 
a nonzero edge for "C eats H" but not for "C is eaten 
by H" in their treatment of a trophic networks). Also, 
note that this linearized system has infinite time con- 
stants at the nodes, i.e., zero values along the diagonal. 
Thus these early models do not include the finite time 
constants that we argue are so important to system dy- 
namics. Later work remedied this omission; the early 
models such as Eq. ^ did not include terms that re- 
searchers subsequently found to be essential for modeling 
real biological systems, such as saturation effects arising 
from resource limitations |14j . Including these additional 
terms leads to a 2 x 2 system matrix that is fully pop- 
ulated with (generically) nonzero terms on and off the 
diagonal. This implies that the resulting linearization 
features finite time constants at each of the nodes, and 
the network is fully connected. That is, where structural 
controllability is concerned, taking into account the full 
dynamics of a food web leads inescapably to the con- 
clusion this system should be controllable with a single 
input. 

DISCUSSION 

Recently, it was reported that sparse inhomogenous 
networks require distinct controllers for a large fraction 
of the nodes to attain structural controllability jTUj . We 
argue that these results are a consequence of assuming 
a special structure for the dynamics at each node: each 
node is treated as a pure integrator. In the application of 
the model set forth in [10] to the real networks considered 
therein, each node is assumed to have an infinite time 
constant. In this paper, we show that (1) for generic, 
arbitrary-order nodal dynamics, structural controllabil- 
ity can be achieved with a single time- varying input, and 
(2) that input should be attached to a PDS. 

The property of a system being controllable has two 
significant interpretations in control theory. First, if a 
system is controllable then it is possible to find an in- 
put to transfer any initial state to any final state in finite 
time. Second, if a system is controllable then it is possible 
to apply a control signal consisting of a linear combina- 
tion of the states that changes the dynamics arbitrarily. 
In particular, it is possible to stabilize an unstable sys- 
tem, a necessary design goal in engineering problems. 
Such a control signal is termed state feedback. 

It is important to note what the first definition of con- 
trollability leaves out. For example, unless the final state 
is an equilibrium, the state will not remain there, but will 
move away. In many engineering applications, it is im- 
portant to find an input that will both stabilize a system 
and hold a specified linear combination (or set of linear 
combinations) of states at desired constant values. This 



is referred to as the problem of setpoint tracking, and 
requires that the system be controllable (so that a sta- 
bilizing control input may be found) and that there are 
at least as many independent control inputs as there are 
linear combinations of states to be held at desired set- 
points [2B]. Hence we see that although one input may 
suffice to achieve controllability of an arbitrary number 
of state variables, in fact the number of inputs limits the 
number of setpoints that may be specified. 

The property of controllability is generically present 
in a system, and thus in practice it is more important to 
know not whether a system is controllable, but whether 
it is almost uncontrollable. In the latter case, the control 
input used to drive the state to its desired value, or to 
achieve the desired dynamics, may be excessively large. 
Hence there is a need for tests — such as those based on 
the control Gramian [TS] — to determine what states are 
almost uncontrollable. In practice these are then treated 
as though they were indeed uncontrollable to avoid the 
excessively large inputs required to control them. 

A more subtle problem arises with the second use of 
the controllability property. In practice, it is rarely pos- 
sible to measure all the states of the system required for 
the control signal used to alter the dynamics of the sys- 
tem. Instead, the control signal is based on estimates of 
the states obtained by processing those states (or linear 
combinations of states) that are measurable. A system is 
said to be observable if it is possible to estimate the states 
using only the available outputs [TS] . As is the case with 
controllability, the property of observability is generically 
present, and it is necessary to determine whether states 
are almost unobservable. 

States that are either uncontrollable or unobservable 
do not influence the input -output relation of a system, 
and cannot themselves be influenced by a control input 
signal based on output measurements. Such systems are 
characterized by a pole (an eigenvalue of the matrix A) 
that does not appear in the transfer function due to be- 
ing canceled by a zero of the transfer function having 
the same value. If the system is almost uncontrollable 
or almost unobservable, then the transfer function will 
have a zero very near to a pole. In this case, it is pos- 
sible to design a control signal based on state estimates. 
However, it may be shown using the theory of fundamen- 
tal design limitations [29, 30J that the resulting feedback 
control system will necessarily have a very small stability 
margin, and be sensitive to disturbances and parameter 
variations. Often, the solution to this problem requires 
the introduction of additional control inputs or additional 
measurements. 

In conclusion, the property of controllability, although 
important, is by no means sufficient to assure a well be- 
haved control problem. One might expect this to be 
true since the property is generically present, as is the 
property of observability. The more relevant questions 
are thus whether the system is almost uncontrollable, 
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almost unobservable, or possesses almost pole-zero can- 
cellations. 
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